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Abstract 

The purpose of this paper is to define the concept of multi-Dirac structures and 
to describe their role in the description of classical field theories. We begin by 
outlining a variational principle for field theories, referred to as the Hamilton- 
Pontryagin principle, and we show that the resulting field equations are the 
Euler-Lagrange equations in implicit form. Secondly, we introduce multi-Dirac 
structures as a graded analog of standard Dirac structures, and we show that 
the graph of a multisymplectic form determines a multi-Dirac structure. We 
then discuss the role of multi-Dirac structures in field theory by showing that 
the implicit field equations obtained from the Hamilton-Pontryagin principle 
can be described intrinsically using multi-Dirac structures. Furthermore, we 
show that any multi-Dirac structure naturally gives rise to a multi-Poisson 
bracket. We treat the case of field theories with nonholonomic constraints, 
showing that the integrability of the constraints is equivalent to the integrability 
of the underlying multi-Dirac structure. We finish with a number of illustrative 
examples, including time-dependent mechanics, nonlinear scalar fields and the 
electromagnetic field. 
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1 Introduction 

To put this paper into context, we begin by giving a historical overview of the use 
of Dirac structures and the associated variational principles in classical mechanics. 

Dirac Structures and Lagrange-Dirac Systems in Mechanics. Originally 
developed by Courant and Weinstein [1988]; Courant [1990] and Dorfman [1993], 
Dirac structures unify (pre-)symplectic and Poisson structures and find their inspi- 
ration in Dirac's theory of constraints, which was established by Dirac [1950] for 
constrained systems including constraints due to degeneracy of Lagrangians. 

It was soon realized that Dirac structures play an important role in mechan- 
ics. Indeed, a wide class of implicit Hamiltonian systems can be naturally de- 
scribed in terms of Dirac structures. In this paper, we will refer to these systems 
as Hamilton-Dirac systems. In particular, it was shown that interconnected sys- 
tems, such as LC circuits, and nonholonomic systems can be effectively formulated 
in the context of implicit Hamiltonian systems (Courant [1990]; Dorfman [1993]; 
Courant and Weinstein [1988]; van der Schaft and Maschke [1995]). On the other 
hand, it was shown by Yoshimura and Marsden [2006a] that Dirac structures in- 
duced from distributions on configuration manifolds naturally yield a notion of im- 
plicit Lagrangian systems or Lagrange-Dirac dynamical systems, which allow for the 
description of mechanical systems with degenerate Lagrangians and with nontrivial 
constraint distributions. Finally, as shown by Leok and Ohsawa [2008], there also 
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exists a natural discrete definition of Dirac structures, whicli can be used to derive 
discrete, geometric approximations of the equations of motion for mechanical sys- 
tems. From the viewpoint of mechanical systems with symmetry, reduction theory 
of Dirac dynamical systems was developed by van der Schaft [1998]; Blankenstein 
[2000] and Blankenstein and Ratiu [2004] , while reduction theory of Lagrange-Dirac 
systems was developed by Yoshimura and Marsden [2007, 2009] and Cendra, et al 
[2010]. In this paper, we will define a graded analogue of Dirac structures, referred 
to as multi-Dirac structures, and we will show that they can be used to describe 
the implicit Euler-Lagrange equations of first-order field theories. 

The Hamilton-Pontryagin Variational Principle. Let L{q, v) be a Lagrangian. 
The Hamilton-Pontryagin principle, originally coined by Livens [1919], is a varia- 
tional principle in which both the position coordinates q and the velocity coordinates 
V are varied independently. The relation q = v is then imposed by means of a La- 
grange multiplier p, leading to an action functional of the form 

Siq, v,p) = (L{q{t),v{t)) + {p{t),q{t) - v{t)))dt. (1.1) 

By taking arbitrary variations with respect to q, v and p, we obtain the Euler- 
Lagrange equations in implicit form: 

dL . dL , , 

Pi = Trj, Pi = Trj, and q=v. 
ov^ oq^ 

In the context of variational principles, it was shown by Yoshimura and Marsden 
[2006b] that the Hamilton-Pontryagin principle and the associated implicit equa- 
tions of motion can be described in an intrinsic way by means of Lagrange-Dirac 
structures. Similar variational principles as the Hamilton-Pontryagin principle have 
been used by Bou-Rabee and Marsden [2009] for the design of accurate variational 
integrators. The Hamilton-Pontryagin principle also appears under a different guise 
in the mechanics of continuous media, where it is related to the Hu- Washizu vari- 
ational principle (see Washizu [1968]). We will extend the Hamilton-Pontryagin 
principle to deal with first-order field theories. 

Stokes-Dirac Structures and Infinite-Dimensional Hamiltonian Systems. 

While there is as of yet no comprehensive theory of Dirac structures for classical field 
theories, an important contribution towards this goal was made by van der Schaft and Maschke 
[2002] , who introduced the concept of Stokes-Dirac structure, an infinite-dimensional 
Dirac structure associated with the exterior derivative and Stokes's theorem. This 
structure can be used for describing Hamilton's equations for a class of field theories 
with boundary energy flow. Recently, it has been clarifled by Vankerschaver, Yoshimura and Marsden 
[2010] that the Stokes-Dirac structure can be obtained through symmetry reduction 
of a canonical Dirac structure that is induced by the canonical Poisson structure on 
the cotangent bundle. In contrast to Stokes-Dirac structures, the multi-Dirac struc- 
tures in this paper are covariant in the sense that time and space are treated on an 
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equal footing. Hence, we expect that a multi-Dirac structure induces a Stokes-Dirac 
structure after choosing a 3+1 decomposition of spacetime; we will return to this 
issue at the end of the paper. 

Multisymplectic Structures. Classical field theories can be described using the 
multi-symplectic formalism, where the jet bundle and its dual of a configuration bun- 
dle play the role of respectively the tangent and cotangent bundle in mechanics. The 
dual of a jet bundle can be equipped with a canonical multisymplectic form, which is 
a higher-degree analogue of the canonical symplectic form on a cotangent bundle. By 
means of the Legendre transformation, this form can be pulled back to the jet bundle. 
These multisymplectic forms provide a natural geometric setting for the Lagrangian 
and Hamiltonian description of field theories. As the literature on multisymplectic 
field theories is very extensive, we can only refer to a subset of it. For fundamen- 
tal aspects, see Goldschmidt and Sternberg [1973]; Kijowski and Tulczyjew [1979]; 
Binz, Sniatycki, and Fischer [1988]; Carihena, Crampin, and Ibort [1991]; Sardanashvily 
[1995]; Cantrijn, Ibort, and de Leon [1996]; Gotay, Isenberg, Marsden and Montgomery 
[1997]; Gotay, Isenberg and Marsden [1999]; Bridges [1997]; Echeverria-Enrfques, Muhoz-Lecanda, and Ror 
[2000] and also the references therein. 

Nonholonomic Mechanics and Field Theories. Much effort has been ded- 
icated to mechanical systems with nonholonomic constraints (see, for instance, 
Vershik and Faddeev [1981]; Bates and Sniatycki [1993]; Bloch [2003]; Cendra, Marsden, and Ratiu 
[2001]). Amongst them, it was shown by Bloch, Krishnaprasad, Marsden, and Murray 
[1996] that the Lagrange-d'Alembert principle plays an essential role in formulat- 
ing the equations of motion from the Lagrangian viewpoint of symmetry and re- 
duction. On the Hamiltonian side, it was shown by van der Schaft and Maschke 
[1994] that constrained Hamiltonian systems can be developed from the viewpoint 
of almost-Poisson structures, which satisfy the Jacobi identity only if the under- 
lying constraints are integrable. This was later extended to the Lagrangian con- 
text by Cantrijn, de Leon, and Martin de Diego [1999]. Further, a notion of im- 
plicit Hamiltonian systems was developed by van der Schaft and Maschke [1995] 
and van der Schaft [1998] in the context of Dirac structures (see also Blankenstein 
[2000]). Nonconservative systems with external forces appearing in servomecha- 
nisms were also illustrated in the context of the constrained Hamiltonian systems 
by Marie [1998]. The equivalence of the Lagrangian and Hamiltonian formalisms for 
nonholonomic mechanical systems was demonstrated by Koon and Marsden [1997, 
1998] together with their intrinsic expressions. 

By analogy with nonholonomic mechanics, nonholonomic classical field theories 
have been developed; for instance, Binz, de Leon, Martin de Diego and Socolescu 
[2002] proposed a generalized version of the d'Alembert principle by using the so- 
called Chetaev condition to represent the bundles of admissible reaction forces as- 
sociated to given nonlinear nonholonomic constraints in the context of the multisym- 
plectic formalism for classical field theories. Vankerschaver, Cantrijn, de Leon and Matin de Diego 
[2005] extended this by relaxing the Chetaev-type conditions to allow for constraint 
forms that need not be determined by the constraint forces themselves. A physical 



1 Introduction 5 



example of a nonholonomic field theory appears in Vankerschaver [2007]. For more 
information about nonholonomic field theories, see also Krupkova and Volny [2006]. 

Just as (almost) Dirac structures and almost-Poisson brackets play an impor- 
tant role in the description of nonholonomic mechanics, so too can nonholonomic 
field theories be described in terms of a multi-Dirac structure induced by the mul- 
tisymplectic form and the constraint distribution. This will be yet another class 
of multi-Dirac structures and as in the case of mechanics we will show that these 
multi-Dirac structures are integrable if and only if the underlying constraints are 
integrable. 

As an interesting aside, we will show that the specification of a nonholonomic 
multi-Dirac structure is not sufficient to completely determine the nonholonomic 
field equations, and that an extra assumption needs to be introduced. This as- 
sumption is trivially true in the case of mechanical systems with nonholonomic 
constraints, and illustrates the difference between classical mechanics and field the- 
ory with nonholonomic constraints. An overview of other conceptual differences can 
be found in Krupkova [2005]; Vankerschaver, Cantrijn, de Leon and Matin de Diego 
[2005]. 

Graded Poisson Structures and Higher-Order Dirac Structures. Not only 
do multi-Dirac structures appear naturally as phase spaces for implicit first-order 
field theories, they also have a number of mathematical properties that are interest- 
ing in their own right. As a special case, multi-Dirac structures include the standard 
Dirac structures of Courant [1990] and the higher-order Dirac structures of Zambon 
[2010]. We will also show that to each fc-form Q, we may associate a multi-Dirac 
structure D, which is (in a sense defined below) the graph of $7. The condition that 
D is integrable is then equivalent to dQ = 0. Furthermore, just as a standard Dirac 
structure induces a Poisson bracket on a restricted class of functions, we will show 
that a multi-Dirac structure endows a subset of the space of forms J7*(M) with a 
graded Poisson bracket, where the graded Jacobi identity is satisfied up to exact 
forms. Similar brackets were defined in the multisymplectic literature (see section 4 
for an overview) and we expect that multi-Dirac structures will be a useful tool 
in studying their properties. We also mention that a similar but different graded 
structure was introduced by Bridges [2006]; Bridges, Hydon, and Lawson [2010] in 
their study of first-order field theories. 

Layout of the Paper. In §2, we briefly review the geometry of jet bundles as basic 
mathematical ingredients and especially describe the geometry of the Pontryagin 
bundle M = J^Y ®Z over a bundle t^xy :Y ^ X. for Lagrange-Dirac field theories. 
In §3, we develop implicit Lagrangian field theories in the context of the Hamilton- 
Pontryagin principle on M, which naturally yield implicit Euler- Lagrangian field 
equations. This is a natural extension of the variational formulation of the Hamilton- 
Pontryagin principle for mechanical systems as in Yoshimura and Marsden [2006b]. 
In §4, we develop multi-Dirac structures on M by introducing a natural pairing 
between multivectors and forms and we show that the graph of a multisymplectic 
form (in some suitable sense) yields a multi-Dirac structure. 
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Further, it is shown in §5 that the imphcit Euler-Lagrange equations for field 
theories can be described using multi-Dirac structures; this leads us to the defi- 
nition of Lagrange-Dirac systems. In §6, we also explore the induced multi-Dirac 
structure from nonholonomic distributions for the case of affine constraints and the 
associated nonholonomic Lagrange-Dirac systems for field theories. Furthermore, 
we develop a constrained version of the Hamilton-Pontryagin principle, referred to 
as the Lagrange-d'Alembert-Pontryagin principle to formulate implicit Lagrange- 
d'Alembert field equations, which are equivalent with the nonholonomic Lagrange- 
Dirac systems for field theories. In §7, we demonstrate the present theory by ex- 
amples of nonlinear Klein-Gordon scalar fields, electromagnetic fields as well as 
time-dependent mechanical systems with affine constraints. Finally, conclusions are 
given in §8 together with some remarks for future works. 

Acknowledgements. We are very grateful to Thomas Bridges, Henrique Bursz- 
tyn, Frans Cantrijn, Marco Castrillon-Lopez, Mark Gotay, Juan-Carlos Marrero, 
David Martin de Diego, Chris Rogers and Marco Zambon who kindly provided sev- 
eral very useful remarks. In particular, we thank Marco Zambon for pointing out a 
flaw in Proposition 6.1 in an earlier version of this paper. 

2 The Geometry of Jet Bundles 

In this section, we provide a quick overview of geometry of jet bundles for the treat- 
ment of Lagrangian field theories. Most of the material in this section is standard, 
and can be found in Saunders [1989]; Gotay, Isenberg, Marsden and Montgomery 
[1997] and the references therein. 

Jet Bundles. Let X be an oriented manifold with volume form rj, which in many 
examples is spacetime, and let ttxy : Y ^ X he a finite-dimensional fiber bundle 
which we call the covariant configuration bundle. The physical fields are sec- 
tions of this bundle, which is the covariant analogue of the configuration space in 
classical mechanics. For future reference, we suppose that the dimension of X is 
n-\-l and that ttxy is a bundle of rank N, so that dimy = n + + 1. Coordinates 
on X are denoted x^, // = 1,2, ...,n -\- 1, and fiber coordinates Y are denoted by 
y^, A = 1, ...,N such that a section : X ^ Y of ttxy lias coordinate representation 

The analogue in classical field theory of the tangent bundle in mechanics is the 
first jet bundle J^Y, which consists of equivalence classes of local sections of vrxy , 
where we say that two local sections (pi, (p2 of Y are equivalent at x G X if their 
Taylor expansions around x agree to the first order. In other words, (pi and (p2 
are equivalent if (pi{x) = (p2{x) and Tx(pi = Tx(p2- It follows that an equivalence 
class [(p] of local sections can be identified with a linear map 7 : TxX — )> TyY such 
that Tttxy 07 = Idr^x- As a consequence, is a fiber bundle over Y, where the 
projection ttjiyy '■ ^ y is defined as follows: let 7 : TxX — > TyY be an element 
of J^Y , then 7rjiyy(7) := y. Coordinates on J^Y are denoted {x^,y^,v^), where 
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the fiber coordinates vj; represent the first-order derivatives of a section. They are 
defined by noting that any 7 G J^y is locally of the form 

The first jet bundle J^Y is an affine bundle over Y, with underlying vector bundle 
the bundle L{TX, VY) of linear maps from the tangent space TX into the vertical 
bundle VY, defined as 

VyY = {ve TyY I Tt:xy{v) = 0}, for y^Y. 

Given any section (f) : X ^ Y oi ttxy, its tangent map T^cj) at x S X is an 
element of JyY, where y = (j){x). Thus, the map x 1-^ T^cp defines a section of J^Y, 
where now is regarded as a bundle over X. This section is denoted j^cp and is 
called the first jet prolongation of (p. In coordinates, j^cp is given by 

x^^^{x>^,y^{xn,d,y''{x^^)), 

where du = d/dx'^ . A section of the bundle J^Y ^ X which is the first jet prolon- 
gation of a section (f) : X ^Y \s said to be holonomic. 

Dual Jet Bundles. Next, we consider the field-theoretic analogue of the cotan- 
gent bundle. We define the dual jet bundle J^Y* to be the vector bundle over Y 
whose fiber at y € 1^ is the set of affine maps from JyY to A^+^X, where K^'^'^X 
denotes the bundle of (n -|- l)-forms on X. Note that since the space of affine maps 
from an affine space into a vector space forms a vector space, is a vector bun- 

dle despite the fact that J^Y is only an affine bundle. A smooth section of J^Y* 
is therefore an affine bundle map of J^Y to A""'"^X. Any affine map from JyY to 
A"+-^X can locally be written as 

where d^^^x := dx^ A dx"^ A • • • A dx^^^ is a coordinate representation of the volume 
form T], so that coordinates on J^Y* are given by {x^ , y^ , p'^, p) . 

Throughout this paper we will employ another useful description of J^Y*. Con- 
sider again the bundle A^+^y of (n + l)-forms on Y and let Z C A^+^y be the 
subbundle whose fiber over y € y is given by 

Zy = {ze A^+^y I i„iy,z = for aU v,w e VyY}, 

where denotes left interior multiplication by u. In other words, the elements 
of Z vanish when contracted with two or more vertical vectors. The bundle Z 
is canonically isomorphic as a vector bundle over Y to J^Y*; this can be easily 
understood from the fact that the elements of Z can locally be written as 

z = pd^'^^x + p'Xdy'^ A d"-Xi„ 
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where d^x^ := 9^ j cT^ x. 

From now on, we will silently identify J^Y* with Z. The duality pairing between 
an element 7 G J^^Y and z G Zy^ can then be written as 

In coordinates, (7,2;) = {p^v^ + p)(F'^^x. 

Canonical Multisymplectic Forms. Analogous to the canonical symplectic 
forms on a cotangent bundle, there are canonical forms on J^K*. Let us first define 
the canonical (n + l)-form on J^Y* = Z hy 

@{z){ui,...,Un+l) = z{T7TYjiY*{ui),...,TTTYjlY*{Un+l)) 

= (^y,jiy*^)('"i'-"''"n+i), 

where we have interpreted z £ Z = J^Y* as before as an (n + l)-form on Y, and 
ui, ...,Un+i € TzZ. The canonical multisymplectic (n + 2)-form on J^Y* is 
now defined as 

n = -dG. 

Denoting again (Px^ := j one has the following coordinate expression 

for G: 

G = pf^dy^ A + p 

while il. is locally given by 

n = dy^ A dp'X A -dp A (2.2) 

It is easy to show (see Cantrijn, Ibort, and de Leon [1999]) that Q is non-degenerate 
in the sense that i;t'^ = implies that X = 0. Moreover, is trivially closed: 
dil = 0. Forms which are both closed and non-degenerate are referred to as mul- 
tisymplectic forms. If the degree of is two, is just a symplectic form, while 
if is of maximal degree, i7 is a volume form. More examples of multisymplectic 
manifolds can be found in Cantrijn, Ibort, and de Leon [1999]. 

A form which is non-degenerate but not necessarily closed is also referred to as an 
almost multisymplectic form. On the other hand, in the context of field theory, we 
will often refer to a closed but possibly degenerate form as a pre-multisymplectic 
form. 

Lagrangian Densities and the Covariant Legendre Transformation. A La- 
grangian density is a smooth map C = Lrj : J^Y A^+^X. In local coordinates, 
we may write 

C{^) = L{x>^,y^,v^)d^+^x, 
where L is a function on J^Y to which we also refer as the Lagrangian. 
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The corresponding covariant Legendre transformation for a given Lagrangian 
density C : J^Y — > A^+^X is a fiber preserving map ¥C : J^Y — )• J^Y* over Y, 
which is given by the first order vertical Taylor approximation to C: 



(F/:(7),7') = A7) + ^ 



A7 + e(7-y)), 



e=0 

where 7,7' S J^Y. In coordinates, ¥ C{x^^ , , v^) = {xf^ , , p'^, p) , where 

Multivector Fields. A useful tool in the description of classical field theory is 
provided by multivector fields. The general theory of multivector fields can be 
found in Tulczyjew [1974], Marie [1997] and the references therein, while applica- 
tions of multivector calculus to classical field theory can be found in (for instance) 
Echevern'a-Enriquez, Muiioz-Lecanda, and Roman-Roy [2002] and Forger, Paufier, and Romer 
[2005]. 

For the sake of generality, we will introduce multivector fields first on an arbitrary 
manifold P . Later on, we will specialize to the case where P is a jet bundle, its 
dual or a combination of both. We denote by T'^p the r-fold tangent bundle, that 
is, the r-fold exterior power A'^'(TP) of TP with itself. A multivector field of 
degree r on a manifold P is a section Xj. of the r-fold tangent bundle T'^P. We 
say that is decomposable when there exist r vector fields Xi, . . . , on P such 
that Xt = Xi A • • • A Xr- An integral manifold of an r-multivector field X^ is 
an embedded submanifold 5 ^ P of dimension r such that X^. spans T"^ S at every 
point. Note that the existence of integral manifolds is not automatically guaranteed 
but depends on certain integrability conditions, which we now describe. 

The integral manifolds of a decomposable r-multivector field = Xi A • • • A X^ 
can be described in an alternative way by noting that X^. induces a distribution /S.Xr 
of rank r, which is locally spanned by the r vector fields Xi, . . . ,Xr. The integral 
manifolds of Xr then coincide with the integral manifolds of A^^, so that in this 
case a necessary and sufficient condition for the existence of integral manifolds is 
given by the integrability of A^^. 

Assume now that P is the total space of a fiber bundle ttxp : P ^ X, where X is 
equipped with a volume form 1] and dimX = n + 1. As mentioned previously, most 
of the times P will be a jet bundle, its dual, or a combination of both. For fiber 
bundles, more can be said about multivector fields and their integral manifolds. 
We will often consider (n + l)-multivector fields on the total space P which are 
decomposable and satisfy the following normalization condition: 

It can be shown that the integral manifolds of such a multivector field Xn+i are 
given as the image of a section ip : X P. This is clear from the coordinate 
representation of Xn+i, as we show. Let X have coordinates {x'^), n = 1, . . . , n -|- 1, 
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and choose a system of bundle coordinates {x^, u^), A = 1, . . . , N on P. Then, the 
multivector field Xn+i can locally be written as 



where the coefficient functions C^{x,u) are local functions on P. Now, we consider 
a section (f> : X ^ P with local coordinate representation (x/^,0^(x)). The 

section (j) determines an integral manifold of -^n+i if it satisfies the following system 
of PDEs: 

3 The Hamilton-Pontryagin Principle for Field Theo- 
ries 

In this section, we introduce the Hamilton-Pontryagin action principle for classical 
field theories. Let us first introduce the Pontryagin bundle 

M := J^Y Xy Z, 

as the fibered product over Y of the jet bundle J^Y and its dual Z. This bundle 
plays a similar role as the Pontryagin bundle TQ®T*Q over a configuration manifold 
Q for the case of classical mechanics. Note also that M is a bundle over Y and also 
over X. 



The Generalized Energy Density. Let C : J Y A""*" X be a Lagrangian 
density which is possibly degenerate. We define the generalized energy density 

associated to C by the map £ : M ^ A^+^X defined as 

£{-f,z) := E{j,z)r] = {z,j) - £(7) 

for {"fjz) € M. In local coordinates {x^,y^,v^,p'^,p) on M, the generalized energy 
density is represented by 

£ = E{x^,y^,v^,p,p^^)d^+'x 

= (p + P>^-^(x^?/^,0)d"+V 

where E = p-\- p^v^ — L{x^, y^, u^) is called the generalized energy on M. 

Pre-Multisymplectic Forms on M. Recall the definition of the canonical forms 
e and n := -dQ on the dual jet bundle Z = J^F*, and let ttzm : M ^- Z be the 
projection onto the second factor. The forms Q and can be pulled back along 
T^ZM to yield corresponding forms on M: 

&M ■= T^*zM^ and Q.M ■= t^*zm^- 
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Note that cannot be multisymplectic since it has a non-trivial kernel: for all 
V € TJ^Y we have that iv^M = 0. As a result, we will refer to Qm as the canonical 
pre-multisymplectic (n + 2)-form on M. Whenever there is no possibility of 
confusion, we will omit the subscript 'M', and denote both canonical forms on M 
respectively by G and Vl. Finally, for any Lagrangian density C with associated 
energy density we introduce another pre-multisymplectic (n + 2)-form on M by 

= Q.M + d£. 



The Hamilton-Pontryagin Principle. Using the pre-multisymplectic forms Q^m 
and fig, we establish a Hamilton-Pontryagin variational principle for classical field 
theories which is similar to the expression (1.1) for mechanical systems outlined in 
the introduction. 

Definition 3.1. Consider a Lagrangian density C with associated energy density £. 
The Hamilton-Pontryagin action functional is defined as 



S{^)= / ri^M-S), (3.2) 
Jx 

where ip is a section of ttxm ■ M ^ X. In other words, tp can be written as ip = 
(7, z) where j is a section o/vr^^ jiy : J^Y — )• X and z is a section of nxz : Z ^ X. 

We have defined the Hamilton-Pontryagin principle in terms of the canonical 
form @M- At first sight, this looks very different from the corresponding principle 
(1.1) for mechanics. However, by using the definition (2.1), we can rewrite the 
Hamilton-Pontryagin action functional as 

S{^P)= [ ((z,iV>-(^,7)+A7)), 
Jx 

where we have written the section tp as tp = (7, z). 

A vertical variation of a section ^ = (7,2;) of ttxm : M — > X is a one- 
parameter family of diffeomorphisms r]\ : M ^ M such that r/o is the identity, r]\ 
is the identity outside of a compact subset U of M, and rjx preserves the fibration 
TTXM : M — )• X, or in other words ttxm ° Vx = t^xm- As a result, the composition 
r/A o -0 is a one-parameter family of sections of M. 

On the infinitesimal level, an infinitesimal variation oi ip = (7, z) is a ttxm- 
vertical vector field Vm ■ ^niip — )• VM, defined on the image of ^p in M, and where 
VM is the vertical subbundle of TM defined by 

VM = {we TM I T-KxAiiw) = 0} . 

Let us denote by X^(M) the module of the vertical vector fields on M. 

One says that ip = (7, z) is a critical point of the action (3.2) if SS^tp) = for 
all variations of ip, where 

d 



- A. 

~ d\ 



A=0 



/ rx{QM-£), 
=oJu 
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where U is an open subset of X with compact closure. 

Proposition 3.2. A section ip = (7, z) of M is a critical point of the Hamilton- 
Pontryagin action functional if ip satisfies the implicit Euler- Lagrange equa- 
tions 

r{VM^^e) = Q, for any Vm &^^iM). (3.3) 

Proof. Recall that a variation of tp is given hy ipx = rjx o tp where rjx : M ^ M is 
the flow of a vertical vector field Vm- Then, it follows that 



A=0 



U 



x=oJu 



[ ^*(VmJ {^m + d^)) + [ d(r(VA/J i&M - £))) 
Ju Ju 



U JdU 



u 

for all Vm, where we utilized Stokes' theorem and the fact that the vrxAf-vertical 
vector field Vm '■ M — )■ VM is compactly supported inU C X. A standard argument 
then shows that ^ is a critical point of S if and only if (3.3) holds. ■ 

Coordinate Expressions. Employing local coordinates {x^ , , , p, p'^) on M, 
the action functional is denoted by 

Ju 

= j^{p + Pa^) ^"^'^ -{iP + P>t)d^^'^ - L{x^ r^^x] 

= I {p'a (I^ - <) + Mx^y^<)} d-^^x. 

(3.4) 

Proposition 3.3. The Hamilton-Pontryagin principle induces the local implicit 
Euler- Lagrange equations for equation (3.3), which are given in coordinates by 

dxi^ dxt^ dy^' dv^' ^ ' 

Proof. By direct computations, the stationarity condition for the action S{'ip) is 
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given by 

A 



{« (1^ - <) . . II) . (-IS - H) -^'^ 
^ X {« (1^ - <) ^ {-"^ ^ H) - {-'S - H) ^'^l ---^ 



JdU 



U 



= 0. 

By definition, the variation of vanishes at the boundary of U, namely, Sy^\gu = 0, 
and so it follows that 

for all variations Sy^, 6v^ and 6p'^. Thus we obtain equation (3.5). ■ 

Generalized Energy Constraint. In addition to the implicit Euler-Lagrange 
equations given in equation (3.5), by imposing the generalized energy constraint 

we can naturally recover the covariant Legendre transformation: 

^-^(^'•■'>^)-||"?- ^^-I|- 

Theorem 3.4. The following statements for a section ip : X ^ M of ttxm are 
equivalent: 

(1) ip is a critical point of the Hamilton-Pontryagin action functional (3.2). 

(2) ip*{XM-i ^s) = for all n x M -vertical vector fields Vm on M . 

(3) ip satisfies the implicit Euler-Lagrange equations (3.5) together with the co- 
variant Legendre transformation (2.3). 

4 Multi-Dirac Structures 

In this section, we introduce the concept of multi-Dirac structures as a natural 
extension to the case of field theories of the concept of Dirac structures on man- 
ifolds developed by Courant [1990]. As a motivating example, we show that the 
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graph of an (almost-)multisymplectic form on a manifold defines an almost multi- 
Dirac structure, and we relate the closedness of this form to the integrability of the 
corresponding multi-Dirac structure. Later on, we will return to multi-Dirac struc- 
tures with nonholonomic constraint distributions and show that they are integrable 
whenever the underlying constraint distribution is integrable. 

Throughout this section, M will be an arbitrary manifold. We will also fix 
an arbitrary integer n + 1 which we will refer to as the degree of the multi-Dirac 
structure, to be defined below. In the remainder of the paper, we will consider the 
special case where M is the Pontryagin bundle M = J^Y Xy Z. In this case, n + 1 
will be the dimension of X. 

Pairings between Multivectors and Forms. Let M be a manifold and consider 
the spaces T'M of /-multivector fields on M and A^'M of A;-forms on M. For k > I, 
there is a natural pairing between elements of S G A^M and X € T'Af given by 

(S,^') := i;tS G A^'~'(M). (4.1) 

We now introduce the graded Pontryagin bundle of degree r over M as 
follows: 

Pr = T''M e A'^+^-r^^ ^4^2) 

where r = l,...,n-|-l. Note that Pr also depends on the choice of n. 

Using the pairing in equation (4.1), let us define the following antisymmetric and 
symmetric pairings between the elements of Pr and Pg as follows. For {X^ S) G Pr 
and (^, S) € P^, where r, s = 1, . . . , n -|- 1, we put 

{{{X, S), {X, S))>_ := \ (i;eS - {-iyn;,f) (4.3) 

and 

(((^, S), {X, S)»^ := i (i;eS + (-l)-i;,S) , (4.4) 

each of which takes values in A"'"''^~^'~*(M). Hence, both of these pairing are iden- 
tically zero whenever n -|- 2 < r + s. 

Let Vs be a subbundle of Pg- The r -orthogonal complementary subbundle 
of Vs associated to the antisymmetric paring ((,))_ is the subbundle {Vg)'^''^ of Pr 
defined by 

(y,)^'^ = {(A',S)GP, |(((A',S),(;P,S)»_ = for all {X,i:)£Vs}. 
Note that {Vg)^''' is a subbundle of Pr, and that (V^)"*"'^' = Pr whenever n + 2 < r + s. 

Almost Multi-Dirac Structures on Manifolds. The definition of an almost 
multi-Dirac structure on M mimics the standard definition of Courant [1990] of 
Dirac structures. 
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Definition 4.1. An almost multi-Dirac structure of degree n + 1 on M is a 

sequence of subbundles Di, . . . , Dn+i, where 

Dr C Pr for r = 1, . . . ,n+ I, 

which is {n + 1)-Lagrangian; namely, the sequence of the bundles Dr satisfies the 
maximally (n + 1) -isotropic property 

Dr = (Ds)^'^ (4.5) 
for r, s = 1, n + 1 and where r + s < n + 2. 

When no confusion can arise, we will refer to the sequence Di, . . . , -Dn+i simply 
as D. Later, we will define an integrable multi-Dirac structure as an almost multi- 
Dirac structure that satisfies certain integr ability conditions. 

For the case of classical field theories, it will turn out that only the orthogo- 
nal complements of the multi-Dirac structure of the "lowest" and "highest" order 
(namely, r = 1 and r = n + 1 respectively) play an essential role to formulate the 
field equations. 

Example. Let M be a manifold with a multi-Dirac structure Di of degree 1. The 
isotropy property then becomes 

Di = {D,)^'\ 
where Pi = TM © T*M and the pairing is given by 

(((A',s),(;P,s)»_ = i(i;eS + i;,s). 

for all {X, S), {X, S) € Pi. This is nothing but the definition of an (almost) Dirac 
structure developed by Courant [1990]. Our concept of multi-Dirac structures also 
includes the so-called higher-order Dirac structures of Zambon [2010]. We will return 
to this at the end of this section. 

Multi-Dirac Structures Induced by Differential Forms. Consider an (n+2)- 
form Ojv/ on a manifold M. We will show that the graph of (in the sense defined 
below) defines a multi-Dirac structure of degree n+ 1. This example of a multi-Dirac 
structure will be fundamental in our subsequent treatment of classical field theories, 
where will be the canonical multisymplectic form, but for now can be an 
arbitrary form. Note especially that at this stage ^Im does not need to be closed or 
non-degenerate. 

Proposition 4.2. Let be an arbitrary (n -\- 2)-form on M and define the fol- 
lowing subbundles Di, . . . , Dn+i, where Dr C Pr '- 

A. = {(^,S) eP. I \x^M = ^} (4.6) 

for r = l,...,n-|-l. Then D = Di, . . . , Dn+i is a multi-Dirac structure of degree 
n +1 on M . 
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Proof. To prove that D = Di, . . . , Dn+i is a multi-Dirac structure of degree n + 1, 
we need to check the isotropy property (4.5), namely namely, Dr = D^'^ for all 
r, s = 1, . . . , n + 1, with r + s < n + 1. 

Let us first show that Dr C Z?^'''. Let {X,T.) G Dr and {X ,T.) € Ds- Then, it 
follows that 

= i{i;eS + (-ir+ii;,s} 

= 0, 

since \^\x^m = (— l)''*!^'!^'^*/- Thus, Dr C -D^"'^. 

Next, let us show that D^'^ C Dr- Let (A", S) G D't'^ ■ By definition of D^'^ , 

i;eS + (-ir+ii;,s = o 

for all (A', S) G D^, i.e. G T'"Z such that ia'^a/ = Then, it follows 

= i;,{(-iri;e^^M + (-ir+'s} 

= 

for all X G T^'Z and with r + s < n + 2. Therefore, one has 

Thus, Dt' C A-. Finally, we have shown that Dr = Dg' ■ So, it follows that 
Dr is a multi-Dirac structure of degree r on M. ■ 

Proposition 4.3. Let D = Di, . . . , -Dn+i a multi-Dirac structure of degree n + 1. 
For any {X, S) G Dr and {X, T,) G Ds, the following relation holds: 

i;es - i-iynxt = 0. (4.7) 

Proof. This is clear from the r-Lagrangian (maximally r-isotropic) property of 

Dr. ■ 

The above r-Lagrangian property of the multi-Dirac structure Dr in equation 
(4.7) may be understood as the field-theoretic analogue of the virtual work prin- 
ciple in mechanics and is related to Tellegen's theorem in electric circuits. We 
shall return to this relation in the construction of Lagrange-Dirac field theories. 
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Wedge Products on Sections of Graded Pontryagin Bundles. We can in- 
troduce a wedge product between elements of the space of sections Pr and Ps , which 
is given by, for {X, S) G Pr and € Ps, 

A(;e,s) := (a'a;P,(((^,s),(a^,s)»^) 

' - 1 . -^^ ^^-^^ 

where on the right-hand side we have used the usual wedge product of multi- vector 
fields. 



Multi-Courant Brackets of Forms and Multi- Vector Fields. We now in- 
troduce a notion of multi-Courant brackets which is a natural extension to the case 
of field theories of the bracket used by Courant [1990]. To do this, we define first 
the Lie derivative £x^ of an /-form S with respect to a fc-multivector field X by 
means of Cartan's magic formula as 

£x^ := i;tdS - {-ifdix^. 

Note that £x'^ is an (/ -|- 1 — A:)-form (see Tulczyjew [1974] for more information). 

Using this notion of Lie derivative, we define a multi-Courant bracket on the 
space of sections of Pr x P, as a graded anti-symmetric bracket 

:r(p,)xr(p,)^r(p,+,„i) 

given by 

[(A',S),(^,S)]^^^ 

:= ([x,x], i;;,s-(-i)(-i)(-i)i:^s + i^'d(((A',s),(;P,s)»^) 
= {[x,xi i:;,s-(-i){-i){-i)i:^s + ^'d(i;eS + (-i)^M;,s)) . 

In the above, [X^ X] is the Schouten-Nijenhuis bracket of the r-multivector field X 
and the s-multivector field X, so that [X, X] is an (r + s — l)-multivector field. 

For the case in which n = and r = s = 1, the multi-Courant bracket cor- 
responds to the standard Courant bracket on the space T(P) of sections of the 
Pontryagin bundle P = TM T*M. 

Definition 4.4. An almost multi-Dirac structure D = Di, . . . , Dn+i of degree n+1 
on M is said to be integrable if for all {X,T,) € Dr and {X,T,) € Dg, where 
r, s = 1, ...,n + 1 and r + s < n + 1, the following condition is satisfied: 



[(^,S),(A',S)]^_^GD,+,_i 



4 Multi-Dirac Structures 18 



Integrable Multi-Dirac Structures. In the construction of the almost multi- 
Dirac structure in Proposition 4.2, we only used the fact that the pre-multisymplectic 
form Qm maps multivectors X into forms S = ix^M- The fact that df^jv/ = was 
left out of consideration, and in fact any (n + 2)-form can be used to define a multi- 
Dirac structure. We will now show that the multi-Dirac structure defined by an 
(n + 2)-form is integrable if and only if is closed. 

Let D = Di, . . . , Dn+i be the multi-Dirac structure defined by an (n-|-2)-form O 
as in (4.6). For D to be integrable, we need to check that the multi-Courant bracket 
(4.9) is closed. Note that the wedge product is always closed: for all {X, S) € 
and {X, T,) € Dg, we have 

(A',s)A(;P,s) = {x,ixnM)Aix,ij^nM) = {XAX/i^ix^M) = iXAX,i;^^j^nM), 

which is an element of D^+s- We now need to check the closedness of the multi- 
Courant bracket (4.9); namely, 

[(^,s),(;P,s)]^^^eD.+._i. 

This condition turns out to be equivalent to the closedness of 1), as is shown in the 
next theorem. 

Theorem 4.5. Let he an arbitrary {n + 2)-form on M which is not necessarily 
closed and let D be the almost multi-Dirac structure defined by the graph ofVtM- 
Then, we have that D is integrable if and only if ATLm = 0. 

Proof. For the almost multi-Dirac structure induced by an [n + 2)-form the 
multi-Courant bracket can be rewritten as follows. Let {X, S) and {X, S) be sections 
of Dj- and Dg. Using the fact that, for the almost multi-Dirac structure, 

S = \x^M and S = i^^m, 
we obtain after some calculations that 

{{X, S), {X, t)}r,s = {[X, X], £xt - + i^d(i;tS + {-lyn^T.)) 

= {[X,X],^;:^^;^^Qm + (-l)"i;^i;edOA/). 

In order for the right-hand side of this expression to be a section of Dr+s, the 
second term has to vanish. This is precisely equivalent to the requirement that 
d^M = 0. ■ 

Multi-Poisson Bracket Induced by a Multi-Dirac Structure. Associated 
to every Dirac structure there is a Poisson structure on a restricted class of func- 
tions, which satisfies the Jacobi identity if and only if the Dirac structure is in- 
tegrable. This correspondence was investigated by Courant [1990] and forms the 
basis for the constrained Dirac brackets of van der Schaft and Maschke [1994] and 
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Cantrijn, de Leon, and Martm de Diego [1999]. We now show that multi-Dirac struc- 
tures similarly induce a notion of multi-Poisson structure, which is a graded ana- 
logue of the usual concept of Poisson structures. Multi-Poisson brackets are studied 
almost exclusively in the context of multisymplectic structures (see among others 
Cantrijn, Ibort, and de Leon [1996]; Castrillon-Lopez and Marsden [2003]; Forger, Paufler, and Romer 
[2005] and the references therein, or Baez, Hoffnung, and Rogers [2010]; Rogers 
[2010b] for a category-theoretic approach) and no comprehensive general theory 
exists as of yet. 

Let D = Di, . . . ,Dn+i be a multi-Dirac structure of degree n + 1. For k = 
0, . . . , n, we say that a /c-form S is admissible if there exists an (n + 1 — k)- 
multivector field Xy, such that 

(A's,dS) 

We denote the space of admissible /c- forms by ^^adm(-^)- define the multi- 

Poisson structure on J^*^^^ associated to D as the map 

{■' '} • ^adm ^ ^adm ~^ ^adm 

(where k,l = 0, . . . , n) given by the following prescription. For S € J^adm ^"^^ 
^ e we put 

{S,S}:=i;,^dS. (4.10) 

where Xj-, is chosen so that (^£,dS) € Dn+i-i- It can easily be shown that the 
bracket does not depend on the choice of multivector Xj^: let X^ be any other 
multivector such that (<Yg,dI]) € The difference {Xj^ — <^g,0) is then also 

an element of Dn~i and therefore we have that 

<((^,S),(^g-4,0)»_=0 

but this is equivalent to 

which shows us that the bracket {S, S} depends only on the choice of E and E. 
At this stage, we remark that the designation "Poisson" is chosen to emphasize the 
similarity with standard Poisson structures. However, to obtain a (graded) Poisson 
algebra we would also have to introduce a product on ^*^^^ with respect to which 
the Poisson bracket is a (graded) derivation. A candidate graded product is given 
by the wedge product (4.8), but a detailed study of these structures is deferred to a 
forthcoming paper. 

To finish this paragraph, we remark that for a multi-Dirac structure induced 
by a multisymplectic form 17, the bracket {E, E} agrees with the multi-Poisson 
brackets of Cantrijn, Ibort, and de Leon [1996]; Forger, Paufler, and Romer [2005]; 
Baez, Hoffnung, and Rogers [2010]. Indeed, let D = Di, . . . ,Dn+i be a multi-Dirac 
structure induced by a multisymplectic form. A A;-form E is admissible if there exists 
an (n -|- 1 — /)-multivector field X-£ so that {X-£, E) € D^+i-i- This is equivalent to 



ixs^Ai = dE. 
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Forms with this property are said to be Hamiltonian. On this class of forms, the 
bracket (4.10) is given by 

{j:,t} = ix^ix^nM, (4.11) 

which agrees (up to sign) with the definition of the brackets in the multisymplectic 
literature. 

It is worth noting that one often encounters two different, non-equivalent multi- 
Poisson brackets in classical field theory. The first one, referred to as the semi- 
bracket in Baez, Hoffnung, and Rogers [2010] is defined as in (4.11), while the 
hemi-bracket is defined for two Hamiltonian forms as {S, S}' := £x^T,. However, 
this notion of bracket does not seem to be well-defined for arbitrary multi-Dirac 
structures, since the right-hand side depends on the choice of multivector field X-^. 
In order words, only the semi-bracket can be defined for multi-Dirac structures. As 
shown by Zambon [2010], a similar phenomenon appears in the theory of higher- 
order Dirac structures, to be discussed below. 

Finally, Baez, Hoffnung, and Rogers [2010] have shown that on a multisymplec- 
tic manifold (M, O) where the degree of Cl is three, the space i^adm '^^^ struc- 
ture of a Lie 2-algebra. This observation was further extended to the case of 
higher-order Dirac structures by Zambon [2010] and to the case of multisymplectic 
structures of arbitrary degree by Rogers [2010a]. Following a conjecture of Zambon 
[2010], it follows that our spaces Q^dmi^^) °f highest-degree admissible forms are 
Lie (n + l)-algebras, but more is true. When considering the space of admissi- 
ble forms ^^adm("^) ™ entirety, we can define a graded Lie bracket on it as in 
Cantrijn, Ibort, and de Leon [1996] as follows. For a € ^^adm(^)' 

|a| = n — k, 

and we define n^^^iM) := ^2dm{M) for A; = 0, . . . , n- 1. The multi-Poisson bracket 
is then an operator 

{., .} : O^dm(M) X ni^^{M) ^ 0^+i(M). 

For multisymplectic manifolds, Cantrijn, Ibort, and de Leon [1996] show that 
is endowed with the structure of a graded Lie algebra, where the graded Jacobi 
identity is satisfied up to closed forms. Their proof can be easily adapted to the 
case of multi-Dirac structures, but we defer a comprehensive discussion of these 
results to a forthcoming paper. 

Relation with Higher-Order Dirac Structures. We now show that multi- 
Dirac structures are naturally related to higher- order Dirac structures as defined by 
Zambon [2010]. A higher-order Dirac structure of order p on a manifold M is an 
isotropic subbundle of TM © A^'(M) equipped with the following pairing: 

{{{X,a),{Y,p))) = ixP + iYa. (4.12) 

and which is involutive under the Courant-Dorfmann bracket: 



{{X, a), {Y, /3)]] = {[X, Y], £x(3 - iya). 



(4.13) 
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Let D = Di, . . . , Dn+i be an integrable multi-Dirac structure of degree n + 1. 
Recall that the multi-Courant bracket is a map from Dr x Ds to Dr+s-i, so that 
Di is closed under the multi-Courant bracket in its own right. We claim that Di 
is a higher-order Dirac structure of degree n + 1. First of all, note that Di d Pi = 
TM © A""'"-'^(M). The antisymmetric pairing (4.3), restricted to Pi, is given by 

(((^,s),(;e,s)»_ = i(i;eS + i;,s), 

which agrees with (4.12) up to an insignificant multiplicative factor of 1/2. Secondly, 
the multi-Courant bracket restricted to sections of Di is given by 

{{X, S) , {X, S)] = {[X, X], - £^T. + d\xT) , 

which is nothing but the Courant-Dorfmann bracket (4.13). 

In this way, we proved that the correspondence between multi-Dirac and higher- 
order Dirac structures, given by mapping to Di is injective. It is 
easy to show that this mapping preserves the Courant-Dorfmann bracket, so that 
integrable multi-Dirac structures are mapped to integrable higher-order Dirac struc- 
tures. Zambon [2010] [Prop. 4.2] shows that this correspondence is a bijection, so 
that multi-Dirac and higher-order Dirac structures are equivalent. 

5 Lagrange-Dirac Field Theories 

We now consider Lagrangian field theories in the context of multi-Dirac structures. 
Recall that X is an oriented manifold with dim X = n + 1 and with a fixed volume 
form rj, locally given by 77 = d^+^x. Recall also that vrxy :Y ^ X \s a. fiber bundle 
and consider a Lagrangian density C on J^Y . Let O be the canonical multisymplec- 
tic structure on Z and recall the pre-multisymplectic (n -|- 2)-form defined on the 
Pontryagin bundle M = J^Y © Z by 17 a/ = t^*zm^- 

Using Qm we can define a multi-Dirac structure D = Di,. . . of degree 

n -|- 1 on M by the construction of proposition 4.2. We refer to D as the canonical 
multi-Dirac structure on M. Explicitly, D = Di, . . . , Dn+i is given by 

Dr = {{X.ix^m) ■■ X € r'M}. (5.1) 

for r = 1, . . . , n -|- 1. 

In this section, we introduce the concept of Lagrange-Dirac field theories, 
which are field theories whose field equations are specified in terms of a multi- 
Dirac structure. More precisely, we shall see that only the component -Dn+i C 
j'n+ij^j X A^(M) of the multi-Dirac structure is needed for the formulation of the 
field equations. We then show that the implicit Euler-Lagrange equations (3.3) for 
an example of a Lagrange-Dirac field theory, as they can be described in terms of the 
canonical multi-Dirac structure D. In the next section, we will address the case of 
field theories with nonholonomic constraints. These field theories can be interpreted 
as Lagrange-Dirac systems as well, where now the nonholonomic constraints are 
incorporated in the specification of the multi-Dirac structure. 
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Partial Multi- vector Fields. Let X be an {n + l)-multivector field on M. We 
say that Af is a partial multivector field if X satisfies 

T^+'ttzm oX = 0, 

where T^~^^ttzm ■ T^'^^M — )• T^+'^Z is the natural projection. In other words, X is 
partial if it does not have a component along the Z-direction. In the remainder of 
this section, we will always consider partial multivector fields which are integrable, 
decomposable and satisfy the normalization condition 

See section 2 for a definition of these properties. Locally, such multivector fields can 
be written as 




where C^, C^^ and are the local component functions of <Y on M. Note that 
the component of the multivector field along d/dv^ is zero. 

Lagrange-Dirac Field Theories. We now define a special class of field theories, 
whose field equations are specified in terms of a Lagrangian £ and a multi-Dirac 
structure D on the Pontryagin bundle M. 

Definition 5.1. Let D = Di, . . . , Dn+i be a multi-Dirac structure of degree n + 1 
on the Pontryagin bundle M and consider a Lagrangian density L with associated 
generalized energy E given by (3.1). A Lagrange-Dirac system for field the- 
ories is defined by a triple {X,E,Dn+i), where X is a normalized, decomposable 
partial vector field, so that 

{X,i-ir+^dE) eDn+i. (5.3) 

Note that in the definition of a Lagrange-Dirac system, only the highest-order 
component Dn+i of the multi-Dirac structure appears. 

At this point, we stress that the multi-Dirac structure D is arbitrary. How- 
ever, when D is the canonical multi-Dirac structure induced by SIm; the resulting 
Lagrange-Dirac equations (5.3) are nothing but the implicit Euler-Lagrange equa- 
tions (3.3) obtained from the Hamilton-Pontryagin variational principle, as we now 
show. 

Theorem 5.2. Let D = Di, . . . , Dn+i be the canonical multi-Dirac structure on M 
induced by VLm o,s in (5.1). The Lagrange-Dirac system (X , E, Dn+i) in equation 
(5.3) induces implicit Euler-Lagrange equations for field theories: 

i^nM = {-ir+^dE, iA^*xMV) = 1, (5.4) 
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which are written in local coordinates as 
together with 

-^{P + P>^-L)=0. (5.6) 
Proof. Recall the pre-multisymplectic form is locally given by 

= dy^ A dp^ A (Tx^j, - dp A 

Note that the fiber coordinates y^, v^,p'^,p appear at most twice in this expression, 
so that the contraction of ^Im with three or more vrxAf-vertical vector fields vanishes. 
Now, the differential of the generalized energy 

Eix^, y^, v^p^p) =p + - L{x^ y^, <) 

on M is given by 

^ ^ dE ^ ^ dE . dE . dE ^ ^ dE ^ 

ox^^ oy^ ov^ ^ op^ op 

-S) + (-1?) + - ^ ^ + *■ 

The local expression for the (n + l)-partial multivector field X is given by (5.2), 
and a computation in local coordinates (see section A. 2) shows that 

xx^M = (-1)"+' [{C^C^x - Ctc\^ - C^)dxf^ + C^dp'X + C^// + dp] . (5.7) 
As a result, we have that (5.4) holds if and only if 

A A AX 

as well as 

Consider now an integral section : X — )■ M of Xn+i, locally represented by 
il){x) = (xf^ , y^ (x) , y^ (x) , p'^{x) , p{x)) . From the equations above, it follows that ip 
satisfies the following system of PDEs: 

^^^A ^^^u 9p _ 
dx'' ^' dx'' ^x^' ^ 

where the local component functions C^, C^^ and are given by (5.8) and (5.9). 
After some simple manipulations, it can then be seen that ip satisfies the following 
equivalent set of equations: 

dp^^dL_ dy^^ J, ^^dL_ 
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together with 



L] 



0. 



These are precisely the equations hsted in the theorem statement. 



The equation (5.6) can be integrated to E = p + p'^v^ — L = a, where a 
is a constant. Without loss of generality we may impose the generalized energy 
constraint by setting a = to obtain 



E = p + p^aVu-L 



0. 



From this and (5.5), we recover the covariant Legendre transformation 



L 



dL ^ 



Pa 



dL 
dvf^ ' 



Matrix Representation. Associated with the implicit Euler-Lagrange equations 
in (5.4), the Lagrange-Dirac system given in equation (5.4) can be also denoted by 

nijX = dE, ix{7r*xzV) = 1, 

where a bundle map Q^j^j : T'^'^^M — )> A^M is defined by the pre-multisymplectic 
structure Q-m by, for every X € T"~^^M, 

ix^M = ^m('^)- 
Then, one obtains a matrix representation as 



Pa 



dL 
_ dL 



V 



\ 



/ 





1 

V 1 



-1 


\ 




\ 























dxi'PA 







/ 


\ dxt^p 


/ 



Note that the above matrix representation for the Lagrange-Dirac formalism is 
related to the matrix representation of the multisymplectic formalism in Bridges 
[1997]; Marsden, Patrick, and Shkoller [1998]. 

Theorem 5.3. The Lagrange-Dirac system {X , E, Dn+i) satisfies the condition of 
the conservation of the generalized energy as 



XldE = 0. 



(5.10) 



In other words, along the solution of the Lagrange-Dirac system {X ^ E, Dn,+i), namely, 
the integral manifold of X, the generalized energy E is constant. 
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Proof. The proof relies on lemma A. 2 for decomposable vector fields. Using the 
field equations, the inner product on the left-hand side of (5.10) can be written as 

XidE = Qm) = 0, 

according to corollary A. 3. ■ 

To see why the previous theorem implies energy conservation, decompose the 
multivector field X as in (5.2). Using lemma A. 2, the interior product can then be 
written as 

k 

XldE = ^(-1)^+^ {X^, dE) 
where is the n-multivector field obtained by deleting X^^ from X, i.e. 

= A 

A=l 

Since the multivector fields X^ are linearly independent, the energy conservation 
equation (5.10) then implies that {X^, dE) = for ^ = 1, . . . , n + 1. In other words, 
the function E is constant on the integral sections of X . 

6 Nonholonomic Lagrange-Dirac Field Theories 

In this section, we develop the idea of nonholonomic Lagrange-Dirac field theo- 
ries. First, we review the theory of mechanical systems in the presence of non- 
holonomic constraints, in order to get an idea of the corresponding structures for 
field theory. Secondly, we define a multi-Dirac structure with nonholonomic con- 
straints by analogy with the induced Dirac structure in nonholonomic mechanics 
as in Yoshimura and Marsden [2006a, b]. We focus especially on the case of affine 
constraints, in which case the induced multi-Dirac structure takes on a particularly 
elegant form. Then, we show how the nonholonomic Lagrange-Dirac field theory 
can be developed in the context of the induced nonholonomic multi-Dirac structure. 

Lagrange-Dirac Systems in Nonholonomic Mechanics. Before going into 
details on nonholonomic Lagrange-Dirac field theories, let us make a brief review on 
nonholonomic mechanics in the context of Dirac structures. 

Let Q be a configuration manifold and let M = TQ ®T*Q be the Pontryagin 
bundle over Q. Recall from Yoshimura and Marsden [2006a] that a Dirac structure 
Daq on Z = T*Q can be induced from a constraint distribution Aq on Q, where 
Aq is not integrable in general, namely, the constraint is nonholonomic. Let L be 
a Lagrangian on TQ, possibly degenerate and let E{q,v,p) = {p,v) — L{q,v) be 
the generalized energy on M. Furthermore, we define a partial vector field in this 
context to be a map X : M TM such that 



Ttttq,m oX = 0, 
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where tttq,m ■ M TQ is the projection onto the first factor. In local coordinates, 
a partial vector field X can be written as X{q,v,p) = {q,v,p,q,0,p). 

We can develop the Lagrange-Dirac dynamical system in the context of the in- 
duced Dirac structure on Pontryagin bundle M (see Yoshimura and Marsden 
[2006b] and Cendra, et al [2010]). To do this, let Am C TM be the distribution 
on M defined by Am = (TvTQAf )~-'^(Aq), where ttqm : M — > Q is the Pontryagin 
bundle projection, which in coordinates is denoted by i^qm{<1^v,p) = q. Since a 
pre-symplectic form Qm is defined by using the projection ttzm : M — >• Z as 

we can define an induced Dirac structure -Da a/ C TM © T*M on M by, for m G M, 

DAiim) = { iwm,/3rn) G T^M X T;^M I Wm G AM{m), 

and I3m - ^^m("i) • Wm G AX/(m) }. 
In coordinates {q,v,p) for m G M, one has 

DAM{q,v,p) = {{{q,v,p), (a, 7,/?)) I <i G A(g),a + p G A°(g), /? = 5,7 = 0}, 

where the constraint set Ag defines a subspace of TQ, which is expressed in a local 
trivialization U C Q by A(g) C M" at each point q £ U. 

A nonholonomic Lagrange-Dirac dynamical system is then specified as 

{X,dE)eDA,„ 

from which we obtain 

dE-ri'M-X e AIj and X G Am- (6.1) 

In local coordinates, the nonholonomic Lagrange-Dirac dynamical system in (6.1) 
is represented as 

p- — eA°{q), q = v, p- — =0, q e A{q). (6.2) 

Nonholonomic Constraints for Field Theories. We now show how multi- 
Dirac structures make an appearance in the theory of classical field theories with 
nonholonomic constraints by analogy with nonholonomic mechanics. For simplicity, 
we restrict ourselves to the case of constraints that are affine in the multi- velocities 
— the general case may be dealt with in a similar way as in Vankerschaver [2005]. 

AfRne Nonholonomic Constraints. Consider a fiber bundle ttxy ■ Y ^ X, 
where X is an oriented manifold with dim X = n + 1 and with a fixed volume form r/, 
locally given by r/ = d^^^x. Let Ay be a distribution on Y. By following Krupkova 
[2005], we say that Ay is weakly horizontal when there exists a distribution W 
contained in the vertical bundle VY such that 



Ay © = TY. 



6 Nonholonomic Lagrange-Dirac Field Theories 



27 



When W is the whole of VY, the distribution Ay is nothing but the horizon- 
tal bundle of an Ehresmann connection on Y. More information on weakly hor- 
izontal distributions and their role in nonholonomic field theory can be found in 
Krupkova and Volny [2006]. In coordinates {x^,y^) on Y, the annihilator Ay is 
spanned by k local one-forms (p" that are given by 

y^) = A'Xix, y)dy^ + yl°(rE, a = 1, k < N, 

where the rank of the matrix A'^{x,y) is k. 

A weakly horizontal distribution Ay defines an affine subbundle C of J^Y as 
follows. We recall that an element 7 = € may be viewed as an injective 
linear map 7 = Tx<j) ■ T^X T^^^.)!^ such that Tttxy 07 = Idx^x- The elements of 
C are then the one-jets 7 taking values in Ay: 

7 € C if Im7 C Ay(y). 

Locally, C can be characterized as follows. Since 7 is given in coordinates by 

we have that 7 takes values in Ay iff 

ij'^ix, y, v) = ^^(x, y) + A'Xix, y)v^ = 0- (6-3) 

In other words, C is locally determined by the vanishing of the k independent affine 
functions tp^. 

From now on, we consider affine constraints of the form (6.3) which come from 
a linear distribution Ay. We say that these constraints are nonholonomic if Ay 
is not integrable. In order to incorporate the nonholonomic constraints into the 
context of multi-Dirac structures on the Pontryagin bundle M = J^Y (B Z over Y, 
we introduce a distribution Am along C as follows: 

Am = (r(7ryM)|c)"^(Ay). 

Note that Am is not a distribution on the entire Pontryagin bundle, but only on 
the affine submanifold C. We denote the annihilator of Am by A^^. Note that A^ 
is locally spanned by the following k linear independent forms 

ip'^ = A'Xidy'^ -v^dx''). 

Another way of defining (^", which also works for nonlinear constraints, is by means 

of the vertical endomorphism S on J^Y. We refer to Vankerschaver, Cantrijn, de Leon and Matin de Diego 
[2005] for further details. 

Nonholonomic Multi-Dirac Structures. Now, we show how a multi-Dirac 
structure on M can be induced from the nonholonomic constraint distribution Am- 
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Proposition 6.1. For r = 1, . . . , n + 1, we define the following subbundles Dj\j^j^r C 
Pr = X A"+2"'^(M), given by 

DA,j,r = {{X, S) G I i;,OM " ^ G A"'"'^'(Am), X G AmA^^''^} . (6.4) 

Then -Daa/ = ^Am,1' • • • ' -^Am,"+i ^■^ ^ multi-Dirac structure of degree n + 1. 
Namely, the sequence of bundles D/s.M,r satisfies the (n + \)-Lagrangian property 

DA,„r = (^A„,.)^'^ (6.5) 
for all r, s = 1, n + l;r + s<n + 2. 

Proof. We have to show that (6.5) holds for all indices r, s with r + s < n + 2. Let 
us first check that 

l?AM,r C (Da„,s)^'^ (6.6) 

Let (A', S) G D^nur and (^, S) G Dam,s- By definition, there exist a G /\^~^'^~^ {A°j^j) 
and a G A"^^^''(A^^) such that 

lA'f^M — S = a and i;pOM — S = q. 

Then, it follows that 

(((^,s),(;P,s)»_ 

= i{i;eS + (-i)-+ii;,s} 

= ^ (i;,5^A/ + a) + i-lY'^'ix ih^M + a)} 
= 0, 

since i^ixCIm = {—^Y^ixix^M, together with i^^a = and i^a. We conclude that 
(6.6) holds. 

Next, let us show that the reverse inclusion 

{Da,,,s)^''' C Z)a„,. (6.7) 
holds. Let G (-Dam,^)^''^- % definition of (^)A„,.)^'^ 

i;eS + (-i)-+ii;tS = o 

for all (;f , S) G -Dam,^) i-^- G Ajv/ A T^-^Z such that 

i;,o;-S = aG A"+'"'(A^./)- 

It follows that 
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has to vanish for ah X G Aj\/ A M. We first choose ^ = 0, so that the 
expression between brackets vanishes. In this case, the condition becomes i^^a = 
for ah a € A"'^^~*(A^) and hence X e Am /\ T'^'^M by lemma A.4. Secondly, for 
arbitrary A" G Am A T^-'^M, we then have the following condition: 

ix {^x^M - S} = 0, 

so that the expression between brackets is an element of /\"~''^~^(A^) by lemma A. 5: 

This concludes the proof of the reverse inclusion (6.7). We conclude that -Daa/ — 
Dam^-i • • • 1 -^AA/.n+i is a multi-Dirac structure of degree n + 1 ■ 

We refer to the multi-Dirac structure -Da^/ = Dahj,i, ■ ■ ■ , D^j^i^n+i in (6.4) as the 
nonholonomic multi-Dirac structure induced by the distribution Am C TM. 
Needless to say that D^^.r is a field-theoretic analogue of the induced Dirac struc- 
ture from a given constraint distribution in mechanics. Furthermore, we note that, as 
in the case of field theories without constraints, only the highest component Dam.^i+i 
plays a role in the formulation of the nonholonomic field equations. Zambon [2010] 
considers an analogue of nonholonomic multi-Dirac structures in the context of 
higher-order Dirac structures. 

Nonholonomic Lagrange-Dirac Field Theories. Consider a nonholonomic 
multi-Dirac structure -Da^/ = D^i,,,!-, ■ ■ ■ ■, Daa/ ,n+i- In the following, only the com- 
ponent -DAM,n+i ^ Pn+i = r"+^M© A^M will be used. Explicitly, this component 
is given by 

DAM,n+i = {{X, S) G Pn+i I i;tw - S G A^^, X G Am A T"M} . 

Now, let C = Lrj : J^Y — > A""'"-^(X) be a Lagrangian density with associated 
generalized energy E = p + p'^v^ — L, and let X he a partial vector field on M. 
Recall that X can locally be written as 

•^=A(i|j+c,f^+a4+'^4)' w-™")-!- 

Definition 6.2. A nonholonomic Lagrange-Dirac system for field theories 
with affine constraints is a quadruple {X,E,Am,Daj^j), where X is a partial 
multivector field, E is a generalized energy, Am is a distribution on M , and Daj^[ = 
-Dam,1) • • • iDAj^j,n+i is a nonholonomic multi-Dirac structure, which satisfies 

iX,{-ir+^dE) £ DA,,,n+l. (6.9) 



and 



X G A''^ A. 



(6.10) 



6 Nonholonomic Lagrange-Dirac Field Theories 



30 



Note that the condition (6.10) in this definition does not appear in the case 
of mechanical systems with nonholonomic constraints. From (6.9) we have that 
X £ A A T^M , and for mechanical systems, where n = 0, this implies that S A, 
i.e. (6.10) holds automatically. For general field theories this is no longer the case, 
and (6.10) arises as an independent condition. 

Theorem 6.3. The nonholonomic Lagrange-Dirac system {X, E, -DAji/,n+i) in- equa- 
tion (6.9) induces nonholonomic Lagrange-Dirac equations for field theo- 
ries: 

^gA^'^'A, Ia-Om - (-l)"+'di? G A^, i;t(vrW?) = l, (6.11) 
which can he written in local coordinates as 

and 



g^ip + p'AV,-L) = 0, (6.13) 
together with nonholonomic affine constraints 

where Aq, are Lagrange multipliers. It follows from equation (6.13) that imposing 
the generalized energy constraint E = reads 

P = L- p'2v^ . 

Proof. By the direct computations using local coordinates, we can easily check that 
equation (6.11) yields (6.12) as well as 

dy^dp^ _ dy^dp^ _ dp_ ^ _dL_ _ ^ ^^^^ 
dxt^ dx" dx" dxf^ dxi' Ox^ " ^^^^ ' 

which is equation (6.13). ■ 



Integrability of Nonholonomic Constraints. We shall investigate the inte- 
grability of the nonholonomic multi-Dirac structure. Recall that the constraint 
constraints in this paper are given by an affine subbundle C of J^Y which induces 
the constraint distribution Aj\/ on M and also that the constraints are said to be 
nonholonomic if A^f is not integrable. 

Consider the nonholonomic multi-Dirac structure I^Aa/ = Dam,1i ■ ■ ■ ;^AM,n+i 
of degree n + 1 defined previously, and endowed with the multi-Courant bracket 
given in (4.9). Then, the following proposition shows that integrability is equivalent 
to the space of section of -Da^ being closed under the multi-Courant bracket. 

Proposition 6.4. The nonholonomic multi-Dirac structure -Daa/ = -Daj,j,i, . . . , D^j^^^n+i 
is integrable if and only if the following conditions are satisfied: 



7 Lagrange-d'Alembert-Pontryagin Principle for Field Theories 



31 



• The distribution Am is integrable; namely, the condition G A^/ holds 
forX,XeAM. 

• The pre-multisymplectic structure Qm satisfies ix/\YAzdi^M = for any three 
X,Y,ZeAM. 

Proof. Recall the multi-Courant bracket on the space of sections of Pr>^Ps, namely, 

[•,•1,^, :r(p,)xr(p,)^r(p,+,„i) 

is given by 

For {X, S) G DAj^f,r and S) G DAj^f,s, one has 

KA", s), , = ([A", X], £xf. - + ^d(i;,s + {-lyn^T.)) 

= {[X,X],\x^x\^M + {-lyix^x^^M)- 
It follows that the multi-Courant bracket is closed iff 

[X, X] E Am a T^+^-^M and dO,/ G A"+'"^''"'^(A1,./). 

This is equivalent to the conditions in the statement of the proposition. ■ 

7 Lagrange-d'Alembert-Pontryagin Principle for Field 
Theories 

Let us show how nonholonomic Lagrange- Dirac field equations can be formulated by 
the Lagrange-d'Alembert-Pontryagin principle, which is an extension of the Hamilton- 
Pontryagin principle to the case in which nonholonomic constraints are given. Before 
going into this, let us review Lagrange-d'Alembert-Pontryagin principle in mechan- 
ics. 

Lagrange-d'Alembert-Pontryagin Principle in Mechanics. We consider again 
the case of a mechanical system first, formulated on a jet bundle vr : y — )• X where 
X = M. The coordinate on X is denoted t and represents time. Let Om = '^*zm^ t)e 
the one- form on M that is induced from the canonical one- form O on Z. We can also 
derive the Lagrange-d'Alembert-Pontryagin by using the Lagrange-d'Alembert- 
Pontryagin principle (see Yoshimura and Marsden [2006b]), which is given by 
the stationary condition for the action integral ^{m) of a curve m{t) in M as 
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{GA/(m(t)) • m{t) - E{m{t))} dt = 



(7.1) 



for all 6m{t) G AMi'm{t)) and with the constraint x{t) € Aj\/(m(t)) and the endpoint 
fixed conditions T'KQM{5m{a)) = TT:QM{5m{b)) = 0. 
By direct computations, it follows that 



A=o 



^ / {@M{mxit)) ■ mx{t) - E{mx{t))} dt 



{d0A./(m)((5m, m) — {dE{m), 5m)} dt 



(7.2) 



{r2A/(m)(m, (5m) — {dE{m),5m)} dt 

= j (^\i{'m) ■ m - dE{m),5rnJ dt 
= 

for all 5m{t) € AA/(m(t)) and with the constraint rn{t) G Ajv/(m(t)) and the end- 
point fixed conditions T7rQM{Sm{a)) = TirQ^^dmih)) = 0. 

Thus, we obtain the Lagrange d'Alembert Pontryagin equations as 



r^Af (m)'' • m — dE{m) G A^j(m) and m G Ajv/ (m). 
In local coordinates {q,v,p) for m G M, it follows that 

(55(g,^;,p) =5 /'{(p(t),g(t)) -i?((?(t),r;(t),p(t))}dt 

{L{q{t),v{t)) + (p(t), g(t) - v{t))}dt = 



(7.3) 



(7.4) 



for all variations 6q{t) G AQ{q{t)) and with the constraint q{t) G AQ(g(t)) and the 
endpoint fixed conditions 5q{a) = Sq{b) = 0. It follows from equation (7.4) that we 
can obtain the Lagrange-d'Alembert-Pontryagin equations as in (6.2). 



Nonholonomic Field Theories and Bundles of Reaction Forces. Now, re- 
call the geometric setting of the nonholonomic Lagrange-Dirac field theories. Namely, 
let TTxY ■ Y ^ X he a fiber bundle, where X is an oriented manifold with 
dimX = n + 1 and with a fixed volume form r], locally given by = d'^'^^x. Let 
M = J^Y Z be the Pontryagin bundle over Y. Let Ay be a distribution on Y. 
As before, suppose that the annihilator Ay is spanned by k one-forms that are 
given in coordinates {x^,y^) for Y by 

^'^{x^, y^) = Aldy^ + A'f^dx^', a = 1, k<N, 
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and Ay defines an affine subbundle C of J^Y which is given in coordinates {x^, y'^,v^) 
for Jiy by 

C:={ix>^,y^,v^)eJ'Y \ = + y) = 0, 

a = l,...,k; fi = I, ...,n + 1}. 

Further, recall that the distribution Aj\/ on M is given by using ttym : M — )• y as 

Am = (T(^yA/)|c)"HAy), 

namely, 

Am := {X G TM \ (^", X) =0, a = 1, k} , 
where are k one-forms on M given by 

= A'Xidy'^ -v^dxi"), 

which span the annihilator of Am- Moreover, the r-fold exterior power of A^ is 
defined by 

r 

Alj = f\AM, r = l,...,n + l. 

We now introduce a closely related concept, namely the bundle of reaction 
forces F C A"^"-'^(Af). This bundle is spanned by the k{n + 1) (n + l)-forms <I>^ 
along C given by 

= A rx^ 

= A'Xidy'^ - v^^dx") A d'^Xf,. 

An intrinsic way of defining the bundle F, as well as further properties, can be found 
in Vankerschaver and Martin de Diego [2008]. 

Definition 7.1. A vertical variation Vm of a section ip of ttxm ■ M ^ X that takes 
values in Am C TM defined on an open subset U with compact closure is admissible 

V'*(Va./J ^) = for all $ G F. (7.5) 

Lagrange-d'Alembert-Pontryagin Principle for Field Theories. Consider 
a Lagrangian density C and with its associated energy density £. Recall from Defi- 
nition 3.1 that the Hamilton- Pontryagin action functional is given by 

5(7, z)= / riQM-£), (7.6) 
Ju 



where U is an open subset of X with compact closure and ip = (7, z) is a section of 
TTxM : M ^ X defined on U . 
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By varying the action S with respect to admissible variations of "0 = (7,^;), a 
critical point tp = (7, z) of the Hamilton-Pontryagin action functional satisfies 



- A. 

~ dX 



A=0 
X=OJU 



(7.7) 



u 



u 



for all vertical vector fields Vm on M that vanish on the boundary of ^{U) and such 
that 

^*(Va./J = for all $ € F. (7.8) 

Let ip : X ^ M he a local section represented in local coordinates by 'ip{x) = 
(x'^, y^(x), f^(x),p^(x),p(x)). The variation (7.7) of the action integral is then 
locally given as 

L {« - <) - ir"^ - ^) A-'S-^) '^'] - » 

for all variations Sy^, and Sp'^, where the variation of y"^ vanishes at the bound- 
ary of X, namely, 6y'^\gj^ = 0. Imposing the condition that the variation be admis- 
sible then yields 

Aldy"" = 0. 

From this, we obtain the following equations of motion for implicit Lagrange- 
d'Alembert equations: 

^ - — - A ^ - „^ - ^ (7q) 

where the are Lagrange multipliers, and 

together with nonholonomic affine constraints 

which serve to determine the k Lagrange multipliers A". Finally, by imposing the 
generalized energy constraint £' = we may suppose that p = L — p'^v^. We 
summarize this in the following theorem. 
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Theorem 7.2. The following statements concerning a section ip of ttxm : M —?■ X 
are equivalent: 

(1) is a critical point of the action (7.6) under admissible variations (7.5); 

(2) ip satisfies the implicit Lagrange- d'Alemhert equations (7.9); 

(3) for all vertical vector fields Vm on M along Aa/ such that ijj*{VM-i $) = for 
all $ G F, 

r{VM^ ne) = 0. 

8 Examples 

In this section, we shah demonstrate some examples of Lagrange-Dirac field theories 
in the context of multi-Dirac structures. Namely, we show that the two examples of 
scalar nonlinear wave equations and electromagnetism can be described as standard 
Lagrange-Dirac field theories. Furthermore, as an example of the Lagrange-Dirac 
field theories with nonholonomic constraints, we consider time- dependent mechanical 
systems with affine constraints. 

Example: Nonlinear Wave Equations. Let us consider the scalar nonlinear 
wave equation {Klein- Gordon equations) in as discussed in Bridges [1997]. 

The configuration bundle is ttxy : Y X, with X = and Y = X x xM.. 
Local coordinates are given by {x^,x^) for U C X and {x^,x^,(p) for Y, and hence 
{x^,x^,cl),vo,vi) for J^Y and {x^,x^,cl),p,p^,p^) for Z = J^Y*. The nonlinear wave 
equation is given by 

where (p is a section of ttxy and 1/ : 1" — )• M is a nonlinear potential. The Lagrangian 
for this equation is given by 

and the Lagrangian density is of course denoted by 

jO{x^, x^, 4>, vo,vi) = L{x^,x^, (j), vo,vi)dx^ A dx^. 
Then, the generalized energy is given by 

E{x^,x'^,(j),vo,vi,p,p°,p^) = p-l-p'^vo+p^vi - L(x°,x\(/),i;o,vi). 

(a) Hamilton-Pontryagin Principle: The Hamilton-Pontryagin principle is given, in 
coordinates {x^ ,x^ ,(j),vo,vi,p,p^ ,p^) for M = J^Y®Z,hy 

5 I {p + p^pfi+p^p^i- E{x^,x^,(l),vo,vi,p,p°,p^))dtAdx = 0. 
Ju 



L{x ,x ,<j),vo,vi) = 
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with the condition of vanishing of 54> at the boundaries of U . By direct computations, 
it fohows that 

and imposing the generalized energy constraints = 0, one has 

p = L{x^,x^,(t),VQ,vi) -p% -P^vi. (8.2) 

(h) Multi-Dirac Structures: On the other hand, since the canonical pre-multisymplectic 
3- form on M is given by 

Qm = -dpo A d0 A dx^ + dp^ Ad4> /\ c?x° - dp A dx^ A dx""^, 

one can define a multi-Dirac structure Dm = -Di , D2 of degree 2 associated to VIm 
on M by the standard prescription. We omit the definition of Z?i, since as discussed 
previously, only D2 is important for the formulation of the dynamics: 

Dm,2 = {{X, S) G P2 = T^Af © A^Af I ix^M = S} . 

In this context, recall that the Lagrange-Dirac system for the field theories is given 
by a triple (Af, D^j) that satisfies 

where the 2-vector field X is given by 



d ^ d d d 



/x=0 

Hence, we have 



9xM ^ >'d(f) ^ ''dp'' ''dp 



. A;- = -dE. 



The direct computation using coordinates {x" ,x ,(j),vo,vi,p,p" ,p ) for M leads to 
vo = Co, VI = Ci, + Cl = —, P P -jr- 

0(p OVq OVi 

dL^ dL 

Co + CqCi — C'oCq = — , Cl + CqCi — CqCi = 

where 

^ _ dcf) _ dp" _ dp 



Hence, by simple rearrangements, we can obtain implicit Euler-Lagrange equations 

(8.1) for nonlinear waves, together with 

^ t T I I 1 \ 1\ 

d^^d^^^"" ^'t>.vo,vi)-p vo-p Vi). 

Needless to say, imposing the generalized energy constraint = yields equation 

(8.2) . 
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Example: Electromagnetism. The multisymplectic description of electromag- 
netism can be found, among others, in Gotay, Isenberg, Marsden and Montgomery 
[1997]. Here, we highlight the role of the Hamilton- Pontryagin variational principle 
and the associated multi-Dirac structure. The electromagnetic potential A = Afj^dx^^ 
is a section of the bundle Y = T*X of one- forms on spacetime X. For the sake of 
simplicity, we let X be with the Minkowski metric, but curved spacetimes can 
be treated equally well. The bundle T has coordinates {x'^,Afj_) while J^Y has 
coordinates {x^,Afj,,A^^^). The electromagnetic Lagrangian density is given by 

C{A,j'A) = -^dAA*dA, 

where the Hodge star operator * is associated to the Minkowski metric. 

(a) Hamilton- Pontryagin Principle: The Hamilton- Pontryagin action principle is 
given in coordinates by 



S 



u 



P ' dx- ^'^ A^^"^ 



d'^x, 



where U is an open subset of X and F^^, is the anti-symmetrization F^j^^, = A^^^^—A^^^. 
The Hamilton-Pontryagin equations follow as before: 

P"'" = ^-. ■4.. = §^. ^ = 0. (8.3) 

(b) Multi-Dirac Structure: The canonical multisymplectic form for electromagnetism 
is given in coordinates by 



0^/ = dA^ A dp^"'" A d'^x^ -dp A d''+^x. 

Let D = Di, . . . ,-Dn+i be the multi-Dirac structure induced by this form and con- 
sider in particular the component D^+i of the highest degree. In terms of this 
bundle, the implicit Euler-Lagrange equations can be described by 

{x,{-ir+^dEL)eDn+i, 

where A" is a partial vector field of degree n + 1 given locally by 

' d ^ d d ^ d 



,1=1 



^dxf "^dA^ ^ dp'^^'' ''dp 
and El is the generalized energy density given by 



El{x^ A^, =p + + -/^luF^^r 

As before, a coordinate computation then shows that the coefficients of X are given 
by 
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and 

dL 

+ C'^ Ckx - Cu^Cx = = 0, 

where 

It is easy to show that these equations are equivalent to Maxwell's equations (8.3) 
in implicit form, together with the energy constraint (5.6). 

Example: Time-Dependent Mechanical Systems with AfRne Constraints. 

Let TTxY : 1^ — 7- X be a finite dimensional fiber bundle called a covariant configura- 
tion bundle over a oriented manifold X. In the case of a time-dependent mechanical 
system, we can set X = M and Y = M x Q, where Q is an n-dimensional configu- 
ration manifold. Namely, Y is an (n + l)-dimensional differentiable manifold, with 
local coordinates {t,q^), where t & X indicates the time and are fiber coordinates 
of Y. The first jet bundle of local sections of ttxy is given by J^y = M x TQ, which 
is denoted, in local coordinates {t,q^,v^) for 7 G J^Y, by 7 = dt® [d/dt + v^d/dq^). 
Let C : J^Y A^X be a Lagrangian density, which is denoted by 

C{j)=L{t,q\v')dt, 

where L is a Lagrangian, possible degenerate, which is defined on J^Y and 7] = dt 
is the volume form on M. The dual jet bundle Z(^ J^y*) = r*y ^ r*M x T*Q 
is the vector bundle over Y, with local coordinates {t,q^ ,pt,p^)- The canonical one- 
form on Z is given by = pidq"^ + ptdt and the canonical two-form is given by 
= -d9 = dq^ A dpi - dpt A dt. 

Let us consider a time-dependent mechanical system which is constrained by 
affine nonholonomic constraints, where Ay is a nontrivial distribution on y = M x Q, 
which is given by, for each (t, q) G Y, 

^Y{t,q) = {jV(i) € J^Y^t,q) I (^r(t,g),iV(i)) = O} , 

where cj)r are the k one-forms given by 

Lpr = A\{t, q)dq' + B''{t, q)dt, r = 1, k. 

Then, the distribution Ay induces an affine subbundle C of J^y = M x TQ locally 
given by 

C = {{t,q\v') G J'Y I AUt,qy+B^{t,q) = 0}. 

Let M = Jiy e Z = (M X TQ) © {T*R x T*Q) be the Pontryagin bundle over 
y = Mx(5— 7'X = ]R and let ttym : M — > y be a canonical projection given by 
{t,q'^,v'^,Pt,Pi) I-)- (t, (?*) and define the distribution on M by 

Am = {T7TYMr\nY{C) C TM (8.4) 
and its annihilator A^j is spanned by k horizontal one-forms on M 



^r = Al{t,q){dq' -v'dt), r = l,...,k. 
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Let P = TM © A^M be the graded Pontryagin bundle and {{{X, S), {X, 
in (4.3) exactly corresponds to 

(((^,s),(<^^,s)»_ = i(i;eS + i;,s) 

for (A', S) G P and S) G P. 

Let TT^M : M ^ Z and the induced (multi-)Dirac structure Dm ■= Dm,i C 
TM © T*M on M induced from the affine constraint distribution Am in (8.4) for 
time-dependent mechanics may be given for the case in which n = and r = s = 1. 
This is given by 

Dm = {{X, S) e TM © T*M \ S - ix^m G A^.^ for X G Am} , (8.5) 

where i^M '■= '^*zm^- 

Now the partial vector field X : M ^ TZ is denoted in local coordinates 
{t,q\v\pi,pt) on M = jiy©Z as 

where rj = dt is the volume form on X = M. Then, the vertical lift of the partial 
vector field X : M TM is given hj X = TirxMiX), where Tttxm '■ TM 
TZ is the bundle map associated to the canonical projection ttxm : M ^ Z. In 
coordinates, one has 

^ = I + ^^1^ +^^|- +^^^1;' i-(-XM.) = 1, (8.6) 

where we note that the horizontal components of f* are zero. 

On the other hand, the generalized energy : M — )■ A^X = M is given by 

E{t,q\v\pt,Pi) =pt+py - L{t,q\v') 

and the differential of E is locally given by 

(8.7) 

Thus, the nonholonomic Lagrange- Dirac system for time- dependent mechanics 
is given by a triple {X^ E, D) which satisfies the condition 

{X,dE) G 

which yields the intrinsic nonholonomic Lagrange-Dirac systems as 

ix^ -dEe Al,. 
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In coordinates, one obtains 

■ .j j . dL „ dL 

t = l, q=v, pi--g- = XrA, Pi- — = 0, 

together with the affine constraints 

A:{t,qy + B\t,q)=0 

and the energy conservation 

^^{pt+P^v'-L{t,q\v^))=0. 
Imposing the generahzed energy constraint E = 0, one can recover 

pt = L -piv\ 

9 Conclusions and Future Work 

In this paper, we have developed a notion of multi-Dirac structure on the Pontryagin 
bundle and its associated Lagrange-Dirac field theories, where we allow a given La- 
grangian to be degenerate. We have shown that the Hamilton- Pontryagin variational 
principle for field theories induces implicit Euler-Lagrange equations for fields, which 
is consistent with the standard case of the Lagrange-Dirac field theory. The ideas are 
shown to be a natural extension of Dirac structures and Lagrange-Dirac dynamical 
systems in mechanics. We have further developed the induced multi-Dirac struc- 
ture from distributions and with the associated nonholonomic Lagrange-Dirac field 
theories. Additionally, we have also shown that the implicit Lagrange-d'Alembert 
field equations, which are equivalent with the nonholonomic Lagrange-Dirac sys- 
tems, can be formulated by the Lagrange-d'Alembert-Pontryagin principle. Finally, 
we have demonstrated the Lagrange-Dirac field theories by illustrative examples 
of scalar nonlinear waves (Klein-Gordon equations), electromagnetism and time- 
depedent mechanical systems with affine constraints. 

Poisson Brackets for Field Theories. In the literature on classical field theo- 
ries, a number of different Poisson brackets have been proposed. We have already 
discussed the multi- Poisson bracket of Cantrijn, Ibort, and de Leon [1996] , but other 

noteworthy developments include the spacetime covariant bracket of Marsden, Montgomery, Morrison, and 
[1986] and the work of Forger and Romero [2005] on relating the multisymplectic 
framework to the Peierls-de Witt bracket. We hope to shed further light on these 
various brackets using multi-Dirac structures. 

Discretizations of Multi-Dirac Structures and Field Theories. Leok and Ohsawa 
[2008] recently proposed a concept of discrete Dirac structure in mechanics. At 
the same time, it has been shown by Bou-Rabee and Marsden [2009] that discrete 
versions of the Hamilton-Pontryagin principle can be used to derive very accurate 
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geometric integrators for mechanical systems. We propose to explore the field- 
theoretic version of these discretizations: since multi-Dirac structures and varia- 
tional principles are a natural way to describe field theories even for degenerate 
Lagrangians, this approach is especially promising. Moreover, in related work (see 
Castrillon-Lopez et al [2010]) we show that the field-theoretic Hamilton- Pontryagin 
principle incorporates the adjoint formalism for covariant field theories. By con- 
structing discrete counterparts of these results, we hope to come up with a discrete 
representation of field theories which respects the gauge freedom of the underlying 
theory. 

Space-Time Decomposition of Multi-Dirac Structures. The multi-Dirac 
structures introduced in this paper are covariant in the sense that no distinction is 
made between time and the spatial variables on the base space. In fact, as we have 
pointed out before, the configuration bundle ttxy Y ^ X can be arbitrary and in 
particular X does not have to represent spacetime. 

While this approach has various advantages (see Goldschmidt and Sternberg 
[1973]; Kijowski and Tulczyjew [1979]; Bridges [1997]; Gotay, Isenberg, Marsden and Montgomery 
[1997] for a discussion), not in the least that the resulting structures are all finite- 
dimensional, it is often necessary to single out one coordinate direction as time, 
for instance when considering the initial-value problem for field theories (see Wald 
[1984]). This procedure is referred to as choosing a (bundle) slicing. A compre- 
hensive overview of the theory of slicings and the corresponding instantaneous for- 
mulations of field theories can be found in Gotay, Isenberg and Marsden [1999]; 
Binz, Sniatycki, and Fischer [1988]. 

We intend to investigate the behavior of multi-Dirac structures in the presence 
of a bundle slicing. Given the choice of a slicing, we expect that a multi-Dirac 
structure will induce an infinite-dimensional Dirac structure on the instantaneous 
space of fields. This is similar to the results of Gotay, Isenberg and Marsden [1999], 
who show that the multisymplectic form on the dual of a jet bundle induces a 
(weakly) symplectic form on the space of fields once a slicing has been chosen. 
Among other things, we expect that establishing the correspondence between the 
theory of covariant, multi-Dirac structures and their instantaneous counterparts will 
be important in a number of areas, most notably the theory of Poisson brackets for 
field theories. We refer to Marsden, Montgomery, Morrison, and Thompson [1986] 
and Forger and Romero [2005] for two different approaches to the theory of covariant 
Poisson brackets for field theories. 

Symmetry and Reduction of Multi-Dirac Structures. The presence of non- 
trivial symmetry groups is an essential feature of many classical field theories, such 
as electromagnetism or Yang-Mills theories. More complex examples are the theory 
of perfect fiuids and Einstein's theory of relativity. It is often useful to obtain a 
form of the field equations from which the symmetry has been factored out; this 
is for instance what happens in writing Maxwell's equations in terms of the field 
strength F^^, rather than the vector potential A^. From a bundle-theoretic point 
of view, various approaches of symmetry reduction for field theories have been pro- 
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posed (see Castrillon-Lopez and Marsden [2003]; Castrillon Lopez and Ratiu [2003]; 
Gay-Balmaz and Ratiu [2010] and the references therein). We propose to develop a 
reduction theory for multi-Dirac structures along the lines of the theory of symmetry 
reduction for Dirac structures in classical mechanics developed by Blankenstein and Ratiu 
[2004]; Yoshimura and Marsden [2007, 2009]. In this way, we expect to unify and 
extend previous results. 



The Relation with Stokes-Dirac Structures. In the introduction, we men- 
tioned the concept of Stokes-Dirac structures, which are infinite-dimensional struc- 
tures based on the exterior differential and Stokes theorem that satisfy the axioms 
necessary for a Dirac structure. These structures were introduced by van der Schaft and Maschke 
[2002] as a means to describe boundary control for field theories in an intrinsic 
way. In Vankerschaver, Yoshimura and Marsden [2010], it was shown that Stokes- 
Dirac structures appear through symmetry reduction for a certain canonical infinite- 
dimensional Dirac structure. 

Given the relevance of Stokes-Dirac structures for the control of field theories, it 
is important to investigate the link with the multi-Dirac structures introduced in this 
paper. Based on the discussion of slicings in the context of multi-Dirac structures 
and the results of Vankerschaver, Yoshimura and Marsden [2010], we may conclude 
that a Stokes-Dirac structure can be obtained starting from a multi-Dirac structure 
by first choosing a slicing and then performing (instantaneous) symmetry reduction. 
One can now ask whether these two operations can be reversed: do we obtain the 
same Stokes-Dirac structure by first performing (covariant) multi-Dirac reduction 
and then choosing a slicing? In other words, when does the following diagram 
commute? 



3+1 

Canonical Multidirac structure ^ Infinite-Dimensional Dirac structure 



Covariant reduction 



Instantaneous reduction 



Multidirac structure ^ Stokes-Dirac structure 

A Generalities about Multi- Vector Fields 
A.l Definitions 

Let M be a manifold of dimension n. A A;-multivector field (where k < n) on M is 
a section X of the fe-fold exterior power T^M of the tangent bundle. The module 
of all /c-multivector fields is denoted by '^^{M). 

We define the interior product of multivector fields and forms as follows. Let 
X G X^{M) and a G Vt\M), where k<l. The left interior product XA a is then 
the unique (/ — /c)-form such that 



[X] a){X') = a{X ^ X') 
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for all A;)-multivectors X'. Similarly, the right interior product of a multivector 
G X'^(M) and a form /3 S Q"^{M), with > m, is the unique {k — m)-multivector 
Xl/3 with the property that 

(^L/3)J7 = (/3A7)W 

for all {k — m)-forms 7 on M. 

Further information about multivectors and interior products can be found in 
(for instance) Tulczyjew [1974]; Marie [1997]. From these papers, we quote the 
following two results on the relation between the interior product and the wedge 
product of vector fields and forms, respectively. 

Lemma A.l (Tulczyjew [1974]). Let X € X^{M) and a G 9}{M). Then 

{X hX')la = {Xla) h X' + {-l)^ X ^ {X' la) 
for any multivector X' , and 

XA {aA/3) = (A'La)J/3 + (-l)^a A (-^J/3) 

for any form (3. 

When the multi-vector field is decomposable, the right interior product with a 
one-form can easily be computed by means of the following formula, the proof of 
which follows from the previous theorem. 

Lemma A. 2. Let X he a k-multi-vector field which is decomposable in the sense 
that X = A^t=i ^^fj-' "^here X^ (fi = 1, . . . ,k) are vector fields and consider a one- 
form a. Then 

k 

Xia = Y^{-lY+^ {X,,,a)X^, 
where X^ is the {k — 1) -vector field obtained by deleting Xfj_, i.e. 

n+l 
A=l 

Corollary A. 3. Let X be a decomposable k-multivector and a an l-form, where 
I > k. Then Xl{XAa) = 0. 

Proof. Let /3 = XA a. With the notations of the previous lemma, we have that 

k 

but {Xf^, /3) = for all ^ = 1, . . . , k, since {X^, f3) = ix^i^n+i ^^ol = 0. ■ 
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In this paper, we also let X^^, (where ^ ^ v) he the (k — 2)-multivector obtained 
by deleting both and X^: 



n+l 



X^y — ^ Xx- 



A=l 

If ^ < i/, we then have that 

Xx A X^, = {-ir+^6^xX, - {-iy+^6,xX^,. 

If /i > i^, note that Xfj^i, = Xy^^. The above formula can then be applied. 

The following lemmas are used in the proof of proposition 6.1. In both cases, let 
A be a distribution on a manifold M and denote its annihilator by A° C T*M. We 
use the following notations for the exterior products: 

A'^A = A A • • • A A and f\\A°) = A° A- - - A A° , 

k times I times 

for all k,l. Note especially that /\\A°) is not equal to (/\'A)°, since the latter is 
given by 

(A'a)° = a° aa'-I(m). 

Lemma A. 4. Consider an l-vector field X and let k be a fixed integer, k > I. If 
ixa = for all a & /\'' A° , then 

X € AAT^^'^M. (A.l) 

Proof. The proof proceeds by induction on the degree / of X. When X is an 
ordinary vector field, the hypothesis is trivially true. For multivector fields X of 
higher degree, we focus first on the case where X can be written as X A 3^ (the 
general case then follows by linearity). We have that 

= ixa = iy (ixa) , 

for all a € /\^ A°. Since ixct € /\''^^A°, this implies by induction that y G 
A A T^^^M. In turn, this implies that X = X Ay satisfies (A.l). ■ 

Lemma A. 5. Consider a k-form a and let I be a fixed integer, I < k. If \xa = 
for allXeAA T^'^M, then 

a G A° A • • • A A° . 

^ V ' 

k times 

Proof. The proof proceeds similarly by induction on the degree k oi a. ■ 
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A.2 Proof of (5.7) 

The computation of equation (5.7) is somewhat involved and depends on a number 
of coordinate identities, hsted here. First, recah that X is an (77, + l)-muhivector 
field with local expression (8), and that the multi-symplectic form Q is locally given 
by equation (2.2). 

Recall the volume form r/ on X is locally expressed by d^+^x and denote 



The following contractions will be useful: 

X^^d^x, = {-ir+^6^, (A.2) 

and, for < f, 

X^uS (i"xA = (-l)"+^+"(5,Adx^ + &^.xdx''). (A.3) 

Now return to the calculation of equation (5.7) and the contraction of X with 
is now given by 

X} Qm = X} {dy^ A dp^ A d^'x^) - X} {dp A d'^+^x), 

and the two terms on the right-hand side will be calculated separately. 
We begin with the first term. Using lemma A.l, we have 

X} {dy^ A dpf^ A d"x^) = {Xidy'^)A {dp'X A d'^x^) + {-l)''+^dy^ A {XA {dp''^ A d"x^)). 

Both terms can be calculated using lemma A.2 and using (A.2) and (A.3). After 
some rearrangements, the result is that 

XA (dy^ A dp>X A d^x,) = (-l)"+2 [{C^C^^^ - C^C\^)dx'^ + C^dp'X 

For the second term, we use the same techniques to conclude that 

X} {dp A d"+^x) = (-1)'^+^ [dp - C^dx"] . 

Putting the results for both terms together, we obtain (5.7). 
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